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ABSTRACT 



Asteroseismology of stars that exhibit solar-like oscillations are enjoying a growing interest with the wealth of observational results 
obtained with the CoRoT and Kepler missions. In this framework, scaling laws between asteroseismic quantities and stellar parameters 
are becoming essential tools to study a rich variety of stars. However, the physical underlying mechanisms of those scaling laws are 
still poorly known. Our objective is to provide a theoretical basis for the scaling between the frequency of the maximum in the power 
spectrum (v max ) of solar-like oscillations and the cut-off frequency (v c ). Using the SoHO GOLF observations together with theoretical 
considerations, we first confirm that the maximum of the height in oscillation power spectrum is determined by the so-called plateau 
of the damping rates. The physical origin of the plateau can be traced to the destabilizing effect of the Lagrangian perturbation of 
entropy in the upper-most layers which becomes important when the modal period and the local thermal relaxation time-scale are 
comparable. Based on this analysis, we then find a linear relation between v max and v c , with a coefficient that depends on the ratio of 
the Mach number of the exciting turbulence to the third power to the mixing-length parameter. 
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1. Introduction 

Scaling relations between asteroseismic quantities and stellar 
parameters such as stellar mass, radius, effective temperature 
and lumin osity have been observationally derived by several au- 
thors (e.g. 



jsity have been observationally de rived by several au- 
Kieldsen & Beddmdl995trChaplin et alJ2fMl2009t 



IStello et alJl2009air using ground-based data. More recently, the 
space-missions CoRoT and Kepler confirmed those results by 
providing accurate and homogeneous measurements for a large 
sample of stars fro m red giants to main-sequence stars (e.g., 
iMosser et al.l l2.Q10b . S caling relations are essential to study a 
large set of stars (e.g., iKallinger et al.ll2009b IStello et al.ll2009bl) 
for which, in general, little is kn own, to provide a first order es- 
timat e for mass and radius (e.g., iBasu et aT 1l2010 HMosser et al 



2009) 



2010), or to probe the populations of red giants (IMiglio et al 



Scaling laws can also lead to a better understanding of the 
• i-h ' underlying physical mechanisms governing the energetical be- 
, haviour of modes. In particular, it has been conjectured by 
?H ■ iBrown et al.1 (Il99ll) that the frequency of the maximum of the 
. . .' power spectrum (v max ) scales as the cut-off frequency v c because 
the latter corresponds to a typical time-scale of the atmosphere. 
The continuous increase of detected stars with solar- like oscilla- 
tions has then confirmed this relation (e.g., iBedding & Kieldsenl 
120031; IStello et a l. 2009a). However, the underlying physical ori- 
gin of this scaling relation is still poorly understood. Indeed, v max 
is associated with the coupling between turbulent convection and 
oscillations and results from a balance between the damping and 
the driving of the modes. The cut-off frequency is associated 
with the mean surface properties of the star and the sound speed, 
making the origin of the y max - v c relation very intriguing. 
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As a first step toward an understanding, one has to deter- 
mine which of the damping rate or the excitation rate is the main 
res ponsible for the max imum of power in the observed spec- 
tra. IChaplin et al 1 (120081) . using a theoretical approach, pointed 
out that in the solar case v max coincides with the plateau of 
the linewidth variation with frequency. We will confirm this re- 
sult using observations from the GOLF instrument in the solar 
case. However, several issues remain to be addressed: is v max for 
any star directly related with the observed plateau in the mode- 
widths variation with frequency? In case of a positive answer, 
what is the origin of this relation? The first issue is quite diffi- 
cult to answer as it is expected to strongly depend on the model 
used for the description of the pulsation-convection interaction. 
Nevertheless, CoRoT observations begin to answer this issue and 
several stars (HD49933, HD 180420, HD49385, and HD52265) 
sugg est that v max correponds to the plateau of the damping rates 
(see iBenomar et al.l 120091: iBarban et all 120091: iDeheuvels etail 
l2010HBallot et al.ll201 1, for details). The second step consists in 
determining the main physical causes responsible for the plateau 
of the damping rates and its mean frequency (vp). Subsequently, 
one has to determine a general scaling law that relates the fre- 
quency of the plateau of the damping rates to the stellar parame- 
ters. In this paper, we discuss the first issue then we focus on the 
second issue by deriving a theoretical relation between vp and 
v c . If one accepts the positive answer to the first issue, this also 
provides the scaling relation between v max and v c . 

This paper is organised as follows. In Sect. 2 we present 
the observed scaling law obtained from a homogeneous set of 
CoRoT data and show that the maximum mode height in the so- 
lar power spectrum coincides with a marked minimum of the 
mode-width when corrected from mode inertia. We then point 
out, in Sect. 3, that such a minimum is the result of a destabiliz- 
ing effect in the super-adiabatic region. The relation between vp 



2 
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and v c is demonstrated in Sect. 4, and conclusions are provided 
in Sect. 5. 



2. The observed scaling law 

We use the CoRoT seimolo gical field data to ens ure a homo- 
geneous sample: HD 49933 feenomar et all l2Q09h HD18 1420 
fearban et all l2009h . HD49385 dDeheuvels et all I2010I) . We 



also use the results on HD508 90 dBaudin et all 1201 1|) and on 
HD 18 1907 (ICarrier et al.ll2010l) . a red giant, and the Sun. The 
characteristics of these stars are listed in Table[T] as well as the 
way their fundamental parameter is obtained. 

For the Sun, due to the presence of pse udo-modes above the 
cut-off frequency (e.g. iGarcia etal 1 ll998h . observational deter- 
mination of v c is not obvious. Nevertheless, one can infer a the- 
oretical relation for this frequency a> c = c s /2H p cc gj -yT^ oc 
MR- 2 T^ /2 (e.g., iBalmforth & Goughj[l990l) . where c s is the 
sound speed, H p the density scale height, g the gravitational 
field, M the mass, R the radius, and F e g the temperature at the 
photosphere. When scaled to the solar case, this relations be- 
comes 

-1/2 



M 

Me 



' effo 



(1) 



with v c0 = 5.3mHz, and M ,R , T e g Q the solar values of mass, 
radius, and effective temperature respectively. Note that we will 
assume H p = H p — P/pg with P,p respectively denoting pres- 
sure, density. Thi s is a commonly used approximation (e.g., 
IStello et aT]|2009ah that presupposes an isothermal atmosphere, 
which is of sufficient accuracy for our purposes. 

Using the stars listed in Table[T| and their measured v max , the 
relation between v max and v c is displayed in Fig. Q] It relies on 
two kinds of results: direct observations of v max in the spectrum 
of the star on one hand, and on estimates of the mass (M), ra- 
dius (R) and effective temperature (r e ff) of the star on the other 
hand. The latter are derived from photometric or spectroscopic 
observations but can be derived in some cases from stellar mod- 
elling. Here, M and R must be derived from stellar modelling 
and not from scaling laws as the aim of this work is to establish 
such a scaling law. The strict proportionality (the fitted slope is 
1.01 ± 0.02) is clearly seen from this sample spanning from the 
Sun to a luminous red giant (HD50890). This is in agreement 
with the r esults obtained by several a uthors for main-sequence 
stars (e.g., Bedding & Kieldsenll2003l) . as well as red giants (e.g. 
iMosser et al.l2010 ). The issue is now to assess the physical back- 
ground underlying this relation. 



3. Height maximum in the power spectrum 

In this section, we confirm that the maximum of the power spec- 
trum of solar-like oscillations is related to the plateau of the line- 
width by using solar observations from the GOLF instrument 
and we then discuss the physical origin of the depression of the 
damping rates (i.e. the plateau). 



3. 1 . Origin of the maximum of height in the power density 
spectrum 

We consider the height H of a given mode in the power spectrum, 
which is a natural observable. To derive it , let us first define t he 
damping rate of the modes given by (e.g., iDupre t et al. 20091) 



-W 



2u\{ r (R)\ 2 M 



(2) 



Star name 


T eff (K) 


M/M Q 


DID 


Vmax (//HZ) 


Sun 


5780 


1 


1 


3034 




OOJU 


1 9 


1 A 


i son 

loUU 


HD181420 


6580 


1.4 


1.6 


1647 


HD49385 


6095 


1.3 


1.9 


1022 


HD52265 


6115 


1.2 


1.3 


2095 


HD181907 


4760 


1.7 


12.2 


29.4 


HD50890 


4665 


4.5 


31 


14 



Table 1. Stellar characteristics (from the literature - see ref- 
erences in Sect for the stars used in the comparison 
with the p resent results. For H D49933, v ma x and J e ff are 
taken from Benomar et al. I (120091) : M and R are taken from 
iBenomaret all (l2010l). F or HD181420, v max and r eff are taken 
from Barban et al. R2009I) : M and R are provided by M.-J. Goupil 
(priv ate communication). For HD49385, v max and T e g are taken 
from iDeheuvels et all (l2010h : M and R are provided by M.-J. 
Goupil (privat e communica t ion). F or HD 18 1907, v, 
are taken from lCarrier et al.l(l2010| For HD50890. 
and R are taken from Baudin et al. 1 (1201 11) . 
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Fig. 1. Frequency of the maximum of oscillation power for the 
main- sequence and red-giant stars of Table Q] as a function of 
the frequency cut-off. All quantities are normalized to the solar 
values. 



where w is the angular frequency, W is the total work performed 
by the gas during one oscillation cycle, £ is the displacement 
vector, and M is the mode mass 



M 



M 



\&(R)\ 2 



■ dm . 



(3) 



£ r (R) corresponds to the radial displacement at the layer where 
the oscillations are measured, M is the total mass of the star. 

For stochas tically excited modes, the power injected into th e 
modes is (e.g., ISamadi & Goup il 2001; Belkacem et al.ll200a) 



1 

8At 



(C 2 R + Cj).. 



(4) 



where C\ and C\ are the turbulent Reynolds stress and entropy 
contributions, respectively. We then introduce the height of the 
mode pr ofile in the power spectrum, which is an o bservable, as 
(see e.g. lChaplin et al.ll2005trBelkacem et alJl2006h 



H 



2rf-M 



(5) 
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Fig. 2. Normalized l/(20 2 ) (dashed-dots line) , where = r\ M, 
and n = V M (solid line) computed from s olar seismic data 
from the GOLF instrument dBaudin et al.l2005l) . The normalized 
mode height (H) corresponds to the dashed line. All quantities 
are normalized to unity at the frequency of v = y max . The associ- 
ated m ode masses are computed as described in Bel kacem et al.l 
(I2006h . 



However, it is useful to express H in a form that does not 
explicitly depend on the mode mass (Al). To this end, we note 
from Eqs. (|2]i and © that both the excitation P and the damp- 
ing rate rj are inversely proportional to the mode mass. Hence, 
to disentangle the effect of the driving and damping from the ef- 
fect of mode mass, we introduce the quantities Yl - P M and 
= rj M, independent of mode masses. Then, using Eq. (0, the 
expression of the mode height becomes 



H 



n 

202 



(6) 



Figure [2] displays the variations of H with mode frequency 
as well as its two contributions n and 1/0 2 . One can clearly dis- 
tinguish a maximum for H near v ^ 3.2 mHz that corresponds to 
the v max frequency, n remains roughly constant (efficient driv- 
ing regime) except at high frequency beyond y max which corre- 
spond s to the inefficient driving regime (see ISamadi & Goupill 
1200 ll for details). On the other hand l/© 2 shows a sharp maxi- 
mum and its variation clearly dominates over that of n and con- 
trols the variation of H and the apparition of its maximum. We 
conclude that the maximum of H is determined by the minimum 
of ® 2 and corresponds to the plateau of the line widths. In other 
words, the depression (plateau) of the damping rates r\ is respon- 
sible for the pr esence of a maximum in the power spectrum, in 
agreement with lChaplin et al.1 (120081) . 

3.2. Origin of the depression of the damping rates 



iBalmforthl (1 19921) mentioned that the depression of the solar 
damping rates originate from a destabilising effect in the super- 
adiabatic layer. He also stressed that the plateau of the damping 
rates occurs when there is a resonance between the thermal time 
scale and the modal frequency. 

Followin g these ideas, we use the MAD non-adiabatic pul- 
sation code dDuprej 120021) for computing the solar damping 
rates. This code includes a time-dependent convection treat- 
ment dGrigahcene et al.l feOOS) different from that by Balmforth 
(1992). Nevertheless, we reach the same conclusion (see 



AppendixlAlfor details): the responsible for the destabilizing ef- 
fect is the Lagrangian perturbation of entropy (5S) that exhibits 
a rapid variation mainly in the super-adiabatic layer as well as in 
the atmospheric layers (see Appendix lA.2l and Fig. IA.2t . 

To understand the origin of such an oscillation and illustrate 
the occurrence of the resonance, we consider the super-adiabatic 
layers and we examine the case of a highly non-adiabatic so- 
lution (see lPesneillll984 for the case of a purely radiative en- 
velope). We assume that Lagrangian perturbations of radiative 
and convective luminosities are dominated by perturbations of 
entropy (see Eq. (1A.6I I and Eq. iA.lh . This leads to a second- 
order equation for the entropy perturbations 6S (Eq. (IA.I6I 1. 
see Appendix. IA.3l for the derivation). To obtain a more explicit 
solution for 6S , we further employ the dimensional approxima- 
tion ddL/dr - 6L/H p , so that 



d In T \ c v / \ c v 



with A = ft-i&, (7) 



where c v - (dU/dT) p with U the internal energy, J[ and S are 
defined by 



/ L c d In c v 
1 L^dlnT 



Y(4-^))(i + (^-i)y) 



S = <2 



(8) 



where kj - (dlni</dlnT) p , L c , and Lr are the convective and 
radiative luminosity respectively, T the temperature, if/ is defined 
by Eq. dA.8b . and we have defined the ratio Q such as 



Q 



with 



<\nr 2 pc v TH p 



rati 



(9) 



with co = 2nv, v the modal frequency, r a local thermal time- 
scale, T ra d and r conv the radiative and convective thermal time- 
scales, respectively. From Eq. (IA. 12) . the oscillatory part of the 
final solution is (6S/c v ) oc exp|-i j SdlnrJ, which describes 
the oscillatory behaviour of entropy perturbations in the super- 
adiabatic layers. 

As discussed in Appendix B3 (Fig B2 top), all modes in the 
range of interest have a similar negative integrated work, W, at 
the bottom of the superadiabatic layer. This corresponds to a 
large damping at this level in the star. In the superadiabatic lay- 
ers, the entropy's oscillatory behaviour controls the oscillating 
behavior of W. When the pulsation period and thus the wave- 
length of the entropy perturbations are too large (<2 <K 1), the 
destabilizing contribution has not grown enough; the cumulated 
work W increases too slowly. The net result at the surface is a 
large damping. When the period is too small (<2 » 1), the rapid 
oscillation of the entropy perturbation causes a rapid oscillation 
of W which increases and again decreases before reaching the 
surface and the net result at the surface is again a large damp- 
ing. Those two limits correspond to low and high frequencies, 
i.e. to the two branches of 1/0 2 displayed in Fig. [2] A minimum 
damping is then obtained for a period neither too small nor too 
large i.e. Q 1 where the destabilizing contribution nearly but 
not quite compensates the strong damping of the layers below 
the super adiabatic layers. 

The value of Q is illustrated in Fig. [3] for three modes. It 
confirms that the resonance Q 1 occurs in the super-adiabatic 
region for the mode with frequency v v max . Hence from the <3 
definition Eq. (01, one derives the resonance condition 



1 

2nT 



(10) 
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V ~ 1 mHz 

v ~ 3. mHz 

~ 4.4 mHz 



4.0 4.1 

log T 




Fig. 3. Product Q (see Eq. (|9]l) versus the logarithm of the tem- 
perature for three values of the mode frequency, for a solar model 
described in AppendixlAl Vertical dotted lines delimit the limits 
of super-adiabatic gradient (see Fig. IA.2l bottom panel). 



4. Derivation of the scaling law 

We now turn to the relation between the thermal time-scale (t) 
and the cut-off frequency. To this end, we use a grid of stel- 
lar models for masses between M = 1 M and M - 1.4 M 
from the ZAMS to the ascending vertical branch, typical of ob- 
served solar-like pulsators. The grid is obtain e d by using the stel- 
lar ev olution code CESAM2k (lMorellll997l iMorel & Lebretonl 
l2008h . The atmosphere is computed assuming a grey Eddington 
atmosphere. Convection is included according to Bohm-Vitense 
mixing-length (MLT) formalism. The mixing-le ngth parameter 
is a - 1.6. The chemical composition follows lAsplund et alJ 
with an helium mass fraction of 0.2485. All quantities 
are evaluated at the maximum of the super-adibatic gradient, 
which corresponds to the maximum of 6S (see Sect. 13.21 ) and 
the location of the resonance (see Eq. (TTOt). 

From Fig.|4](top), the relation between the thermal frequency 
(1/t) and the cut-off frequency (v c ) is close to linear but still 
shows a significant dispersion. More precisely, the relation be- 
tween those two frequencies is approximatively linear and the 
dispersion is relate d to the dispersion in mass, in agreement with 
observations (e.g. iMosser et al. 1 120101) . We then conclude that 
the observed relation between v max and v c is in fact the result 
of the resonance between y max and 1/t, as well as the relation 
between 1/t and v c . 

To go further, let us investigate the relation between 1/t and 
v c . First, Eq. ^ can be recast as 



pc v TH p 



1 + 



(ID 



where F conv and F m d are the convective and radiative fluxes, re- 
spectively. The MLT solution for the con vective flux and the con- 
vective rms velocity can be written (see ICox & Giulil[l968l for 
details) 




Fig. 4. Top: Thermal frequency (1/t) computed from Eq. (0 
versus the cut-off frequency computed following Eq. ([T), nor- 
malized to the solar values, for models with masses ranging from 
M = 1M Q to M = 1.4M Q (with an increment of 0.05 M ) and 
from the ZAMS to the ascending vertical branch. Bottom: The 
same as for the top panel, except the thermal frequency is di- 
vided by the Mach number to the third versus the normalized 
cut-off frequency. 



where A = aH p is the mixing length, a the mixing-length pa- 
rameter, V = (dlnT/dlnP), V = (dlnT'/dlnP) the gradient 
of rising convective element, 2 = (31np/<91n T)^ P , with // the 
mean molecular weight, and Fi = ((91nP/<91np) ad . Now, by in- 
serting Eq. (TT21 and Eq. (fT3l into Eq. (fTTT i. one obtains 



1 



/ r2 \ 
1 1 



2H, 



1 + 



rad 



(14) 



where M a = v conv /c s the Mach number, and Xp - 
(dinP/SlnpV 

We verified that for a given physic, the ratio F ra d/-Fconv is ap- 
proximately the same for all the models considered in the super- 
adiabatic layer. Hence, by use of Eq. ( fT~4T > as well as the reso- 
nance condition (Eq. (fTUb). we conclude that 



l pc,,v com T-^- (V - V) 



2 

2 V2r| /2 



(V-V) 



1/2 



(12) 
(13) 



1 



oc — oc 
T 



/ r 2 \ 
1 1 



a 



(15) 



which is the observed scaling between v max and v c (see Fig.[TJ, 
since the thermodynamic quantities hardly vary. 
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Equation ( fT5b describes the observed scaling between Vmax 
and v c (see Fig.[T} but also shows that most of the departure from 
the linear relationship between 1 /r and v c comes from the Mach 
number, as confirmed by Fig. |4] (bottom panel). We also point 
out that as shown by Fig. [4] and Fig. Q]for the main-sequence 
stars, the departure from the linear relationship is of the same 
order of magnitude as the uncertainties on the cut-off frequency. 
However, our grid of models is not suited for a proper compar- 
ison between the observations and the theoretical relation. This 
work is definitely desirable in the future. 

5. Conclusion 

We have addressed the issue of the physical reason for the exis- 
tence of a scaling relation between v max and v c . We have found 
that the depression of the damping rates determines v max because 
there is a resonance between the local thermal time-scale in the 
super-adiabatic region and the modal period. This implies that 
^max does not scale only with v c but also with the ratio Ail/a. 
As pointed out in Sect.Q] the observed scaling between v max and 
v c is not obvious at first glance since the first frequency depends 
on the dynamical properties of the convective region while the 
second is a static property of the surface layers. The additional 
dependence the Mach number resolves this paradox. 

This scaling relation is potentially a powerful probe to con- 
straint the dynamical properties of the upper-most layers of 
solar- like pulsators through the ratio Ail /a. Indeed, as shown in 
this paper, most of the dispersion in the v max -v c scaling is related 
to the Mach number. The investigation of the ratio between v max 
and v c in main-sequence stars, subgiants, and red giants may give 
us statistical information on the evolution of the properties of tur- 
bulent convection from main-sequence to red giant stars, through 
for instance the mixing-length parameter. Indeed, a future work 
will consist in computing models, that correspond to the obser- 
vations, and to make a comparison between the observed and 
theoretical dispersion from the linear relation between v max and 

Ve- 
in other specific cases, for which stellar parameters are well 
known (e.g., in pulsating binaries) the relation between v max and 
v c could gives us directly the value of the Mach number in the 
upper-most convective layers. 

Acknowledgements. K. B. gratefully acknowledges support from the CNES 
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Appendix A: The plateau of the damping rates 

A. 1. Computation of the damping rates 

Damping rates have be en computed with the non-adiabatic 
pulsation code MAD dDupretl |2002|) . This code includes 
a time-dependen t conv ection (TDC) treatment described in 
iGrigahcene et al.l (I2005I) . This formulation involves a free pa- 
rameter j6 which takes complex values and enters the perturbed 
energy equation. This parameter was introduced to prevent the 
occurence of non-physical spatial oscillations in the eigenfunc- 
tions. We use here the value ft = -0.55 - 1.7; which is cali- 
brated so that resulting damping rates reproduce the variation 
of the solar damping rates r\ with frequency and more precisely 
the depression of the rj profile (see Fig. IA.U . Note that TDC is 
a local formulation of convection. This simplifies the theoreti- 
cal description and is sufficient here as we seek for a qualitative 
understanding of the relation between the frequency location of 
the damping rate depression and the cut-off frequency. We stress 
that the above approximations do not influence qualitatively the 
conclusions. 

This approach takes into account the role played by the vari- 
ations of the convective flux, the turbulent pressure, and the dis- 
sipation rate of turbulent kinetic energy. Hence, the integral ex- 
pression of the damping can be written as follows 



1 



2o>M\MR)\ 



1 Jo 
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6p' 
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tn rb 



+ (T 3 - l)T6S 



Am 



(A.1) 

where £; r {R) is the radial mode displacement at the photo- 
sphere, o> the mode frequency, p the mean density, r$ — 1 = 
(din r/d In p) s , T the unperturbed temperature and the star 
denotes the complex conjugate. The symbol 6 represents a 
Lagrangean perturbation: 6S is the perturbation of specific en- 
tropy, dp the density perturbation, SP tm t the perturbation of tur- 
bulent pressure. The quantity 6P tuI b/p represents the contribu- 
tion of turbulent pressure while the second term (Tj - 1) T 5S 
includes the variations of radiative and convective fluxes as well 
as the dissipation rate of turbulent kinetic energy, as given by the 
energy conservation equation 



icrT6S 



A5L r A6L C 



Am 



dm 



+ 5e, 



(A.2) 



with 5L r , 5L C being the perturbations of the radiative and con- 
vective fluxes respectively, be, the perturbation of the dissipation 
rate of turbulent kinetic energy into heat, and cr — to + it]. Note 
that Eq. dA.2b is only valid for radial modes we are interested in. 



A.2. Origin of the depression of the damping rates 

The depression of the damping rates, located around v ~ 3.5 
mHz (Fig lA.ll i. results from a subtle balance between the above 
contributions to the work integral. The cumulated work integral 
(regions where it increases outwards drive the oscillation and re- 
gions where it decreases outwards damp the oscillation) allows 
us to identify the processes that create this depression. Fig. IA.2l 
(top) shows that mode damping results from stabilizing effects 
from inner layers at temperature greater than log T - 4 desta- 
bilizing effects in the upper layers located in the super-adiabatic 
layers (i.e. between log T ~ 3.95 and log T ~ 3.8) and for high 
radial order modes again stabilizing effects from the very outer 
layers. Hence, the behavior of the product of the damp- 
ing rates to the mode mass, which is the integral appearing 
in Eq. (1A.1I ). can then be described as follows: for modes with 




Fig. A.l. Product of the damping rates {iff times the mode mass 
(M) versus mode frequency. The vertical dotted lines identify 
radial mode orders. 



frequencies v < v max the higher the mode frequency the larger 
the contribution of the destabilizing region and keeps on de- 
creasing. For v > v max , despite an increasing contribution of the 
superadiabatic boundary layers, atmospheric layers stabilize the 
modes resulting in an increase of 0. At v = v max , compensation 
is maximal giving rise to the minimum of 0. 

The physical cause of the destabilizing effects in the su- 
peradiabatic regions is revealed by Fig. IA.2I (middle). The 
Lagrangian perturbation of entropy exhibits a rapid variation that 
occurs mainly in the super-adiabatic layer and in the atmospheric 
layers. As the frequency of the mode increases, the amplitude of 
this variation (which is a spatial oscillation as seen in the next 
section) also increases. The wavelenth of this spatial oscillation 
decreases with increasing frequency. This causes a similar be- 
havior of the cumulated work. 



A.3. Oscillation of entropy fluctuations 

To understand the behaviour of 6S in this region, let us first ex- 
amine the fluctuations of radiative and convective luminosity ap- 
pearing in the energy equation (Eq. (1A.2» . 

In the diffusion approximation, the fluctuations of radiative 
luminosity reads 



„&^„6T 6k Sp A5T/Ar A% r 

Z h 1 ~ A-r l A ~A \ A -->) 

r Tup AT/Ar Ar 



where £ r is the mode's radial displacement, 6T the Lagrangian 
perturbation of temperature, 5k the perturbation of opacity, 
and k the opacity. By using the perturbed continuity equation 
Eq. ( 1A.3I I. becomes, for radial modes, 



6L R I d I6T\ 6T 6k 
Tr~ ~ AT/Ar A~r \~T) + ~f ~ 
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(A.4) 



where we have neglected f r /r compa red to d$ r /dr. This assu mp- 
tion is valid for radial p modes (see iBelkacem et al.l 120081 for 
details). We further assume that, in the super-adiabatic region, 
perturbation of temperature fluctuations and opacity are domi- 
nated by entropy fluctuations, so that 
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Fig. A.2. Top: Normalized cumulated work integral versus log- 
arithm of temperature for five values of eigenfrequencies. These 
modes are emphasized in Fig. IA.1I Middle: Imaginary part of 
the Lagrangian perturbation of entropy versus logarithm of tem- 
perature. Bottom: Normalized convective and radiative fluxes 
versus logarithm of temperature. The difference between the real 
and adiabatic gradient (V) is also plotted and normalized to its 
maximum. 



where c v - (dU/dT) p with U the internal energy per unit 
mass, and kj = (<91n/e/<91nr) p . Hence, inserting Eqs. ( IA.5t in 
Eq. (IA.3b we obtain 



^ 7 d/*S\ £ 
L R dr/dr dr \ c v ) c v 



(A.6) 



Fig. A.3. Top: Real part of the radiative luminosity perturbation 
for the mode of radial order n — 23 versus logarithm of tempera- 
ture. The solid line represents the full non-adiabatic computation 
as described in Sect. IA.1I while the dashed line represents the 
approximate expression given by Eq. (IA.6t . The rapid variation 
near log T = 4.17 is the result of the presence of a node of the 
eigenfunction. Bottom: Real part of the convective luminosity 
perturbation for the mode of radial order n — 23 versus logarithm 
of temperature. The solid line represents the full non-adiabatic 
computation as described in Sect. IA.ll while the dashed line rep- 
resents the approximate expression given by Eq. (IA.71 . 



The approximate expression Eq. (1A.6K even if imperfect, cap- 
tures the main behaviour of SLr/Lr in the superadiabatic bound- 
ary region, as shown by Fig. IA.3l (top). Note that the disagree- 
ment observed in the inner layers in Fig. IA.3l (top panel) is due 
to the approximation Eq. dA.5t since for those layers the density 
fluctuations are dominant. However, we are mainly interested in 
the super-adiabatic region (log T < 3.9) where Eq. dA.6b is suffi- 
ciently valid for our purpose. 

We now turn to the Lagrangian perturbation of convective 
luminosity. It is dependent on the adopted time-dependent treat- 
ment of convection. Consistent with Sect. IA. II we use the for- 
malism developed by iGrigahcene et al.l (120051) . A good approx- 
imation of their Eq. (18) in the super-adiabatic layer, as shown 
by Fig. lA.3l (bottom). is 
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with 



where 



ilr-C 1 + 
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(A.8) 



where t c is the life-time of convective elements, cjr is the char- 
acteristic cooling frequency of the turbulent eddies, and 
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(A.9) 
(A. 10) 
(A.ll) 



where A = 8/3 is a constant introduced bv lUnnold 19671) to close 
the equation of motion describing convection, and wg is the char- 
acteristic cooling frequ ency of turbulent eddies due to radiative 
losses (see Eq. (C12) of lGrigahcene et al. 20051). 

We are now interested in obtaining the equation that qual- 
itatively explains the oscillation observed in Fig. IA.2I (middle 
panel). Hence, one has to exhibit in an analytical way the fre- 
quency dependence of the entropy fluctuations (6S). To this end, 
we will use two different assumptions. The first, and most im- 
mediate way is to assume that in the energy equation (Eq. dA.2b ) 
ASL/Ar - 6L/H p . This is a crude approximation, but which per- 
mits to immediately exhibit the role of the Q factor. Then using 
Eq. (1A.6I ) and Eq. ( lA.7b . one obtains 



0, with A = $l-i&, (A.12) 
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where J{ and & are defined by 
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where Hj is the temperature scale-height. 

To derive an analytical solution of Eq. (IA. 16b is not triv- 
ial. Hence, further simplifications are needed. We then assume 
the coefficients T, Q, "H are constant. Assuming solutions of the 
form (6S /c v ) cc e klnT , one has the solutions for k 



'1,2 



-§ + [g 2 - mt\ 
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(A. 18) 



At the maximum of the super-adiabatic gradient, the radiative 
luminosity dominates over the convective ones. Hence, we fur- 
ther neglect the ratio L c /L compared with Lr/L. Eq. ( IA. 1 8b then 
simplifies to 
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AG 
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(A.19) 



From Eq. ( 1A.19K one concludes that for Q l,kis real and 6S 
does not oscillate. This corresponds to the limit of low-frequency 
modes for which both k and the imaginary part of SS are small, 
as confirmed by the full numerical computation presented in 
Fig. IA. 21 (middle panel). In contrast, for (J » 1 (i.e., for large 
frequencies) the imaginary part of the wavenumber increases as 
depicted by Fig. lA.2l (middle panel). 

Eventually, both methods to derive the frequency behaviour 
of 6S converge toward the same conclusion, i.e. that the factor <3 
explains the oscillation of entropy fluctuations and its frequency 
dependence. 



where we have defined the ratio <2 as 



Q - lot , with t 1 = 



Anr 2 pc v TH p 



(A. 14) 



with t is a local thermal time-scale. Note that we have neglected 
the imaginary part of <x in Eq. ( IA. 14b - We stress that this thermal 
time-scale can be recast into 



T- 1 = T- 1 + 
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(A. 15) 



where r conv and r la d are associated with the convective and radia- 
tive luminosities, respectively. From Eq. (IA.12K the oscillatory 
part of the final solution is (6S /c,,) oc exp \-i J Sdln rj, which 
explains the oscillatory behaviour of entropy perturbations in the 
super-adiabatic layers and its frequency dependence. 

An alternative way to proceed is to use the energy equation 
(Eq. dA.2b ) together with Eq. ( IA.6b and Eq. iA.li . one obtains the 
second order differential equation 
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